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O 

The use of statistical fluctuations as probes of the microscopic dynamics of hot and 

dense hadronic matter is reviewed. Critical fluctuations near the critical point of QCD 

matter are predicted to enhance fluctuations in pionic observables. Chemical fluctuations, 

especially those of locally conserved quantum numbers, such as electric charge and baryon 

number, can probe the nature of the carriers of these quantum numbers in the dense 

O ■ medium. 
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1. INTRODUCTION 



Thermal models for the production of final-state hadrons in nuclear collisions have been 
extremely successful over a wide range of collision energies ranging from the Bevalac/SIS, 
over the AGS and SPS, to the RHIC The parametric dependences of the temper- 

ature T and baryon chemical potential \i on the nucleon-nucleon center-of-mass energy 
\fs map out two lines in the phase diagram, which are usually called the chemical and 
thermal freeze- out lines. The chemical freeze-out line (T c h, /i c h) is deduced from the abun- 
dance ratios of different species of hadrons, while the thermal freeze-out line (T th ,/i th ) is 
deduced from the slopes of the transverse momentum spectra of the emitted hadrons. 

What do these results tell us about QCD except that the hadronic matter at freeze-out 
is quite well described by a state of thermal and chemical equilibrium? Firstly, rather 
strong arguments have been presented that the chemical and thermal equilibrium can only 
be established by some complex "prehadronic" dynamics, which requires the presence of 
gluons as dynamical degrees of fredom ||. Secondly, the thermal equilibrium description 
requires the assumption of a strong transverse collective flow of the matter at freeze-out, 
which exhibits an increasing degree of azimuthal anisotropy as the CM energy increases. 
This effect, called elliptic flow [f7|, indicates the presence of a strong transverse pressure 
very soon (1 fm/c) after the begin of the nuclear reaction ||. 

As the chemical freeze-out parameters at SPS and RHIC energies are very close to 
the expected phase boundary between hadronic matter and the quark-gluon plasma, it is 
reasonable to conjecture that the produced matter was first thermalized in the deconfined 
phase and then evolved through the phase boundary into a thermal gas of hadrons. The 
crucial question is whether any signatures remain in hadronic observables that carry 
information about this evolution prior to the final freeze-out. 

As I will explain, there are reasons to believe that statistical fluctuations |J in certain 
hadronic observables may preserve information about earlier times, even though their av- 
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erages only tell us about the freeze-out conditions. It is here where the finite volume and 
the finite life-time of the hadronic fireball produced in relativistic heavy-ion collisions turn 
out to be advantages rather than deficiencies. The relative size of fluctuations in ther- 
modynamic observables, as compared to the average values, decreases inversely with the 
volume V of the system. And if the system would evolve infintely slowly, the fluctuations 
at freeze-out would only reflect the freeze-out conditions. 

The talk is organized as follows: After a brief review of statistical fluctuations in parts 
of an equilibrated system, I will discuss fluctuations in the momentum spectra of emitted 
particles, and later fluctuations in particle abundances, or chemical fluctuations, where 
the main emphasis will be on the fluctuations of conserved quantum numbers. The talk 
will conclude with some experimentally relevant considerations and an outlook. 

2. STATISTICAL FLUCTUATIONS IN A SUBVOLUME 

Consider small subsystem of a large system in thermodynamic equilibrium at temper- 
ature T = f3~ l . Let us ask how the average energy (E)y contained in the subvolume V 
changes as a function of temperature. A straightforward calculation yields: 

d _ d TryjHe-f" 1 ) 2 2 

dp { W ~ W Tr v (e-0») ~ ~ {E )V + {E)v - (1) 

As the right-hand side is just the fluctuation of the energy contained in V, we have 
obtained the important result, that the energy fluctuations are determined by the heat 
capacity of the matter in the subvolume: 

(AE 2 ) V = = T*C V (2) 

where Cy is the heat cpaacity. Fluctuations of other (conserved) quantities O can be 
calculated in a similar way, by considering the variation of the grand canonical average 
Ov,t with respect to a change in the chemical potential /i associated with O: 

{A ° 2)v = T ^ff I " TX °- (3) 

The expression on the right-hand side, apart from the factor T, is called the susceptibility 
Xot which measures the response of the medium to a change in the chemical potenial. 
associate with the observable O. An important example is the magnetic susceptibility 
ch%Mi where the observable is given by the magnetization M, and the magnetic field H 
assumes the role of the chemical potential. Examples of interest in the case of hadronic 
matter are the chiral susceptibility Xm and the Polyakov loop susceptibility xl, which 
measure the response of the medium to a change in the quark mass, and to the addition 
of a free heavy quark, respectively. 

3. THERMAL FLUCTUATIONS 



As outlined above, fluctuations of the thermal energy in a given subvolume are a mea- 
sure of the change of the average thermal energy with temperature, i.e. the heat capacity 
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of the matter contained in the subvolume. A similar relation connects the fluctuations of 
the local temperature to the heat capacity Cy\ 

^1 = -L. (4) 
T 2 Cy K ' 

This relation could be used to measure the specific heat of the hadronic matter created in 
high-energy nuclear collisions |10| , |TT|| . The idea here is to use the slope of the transverse 
momentum spectrum of emitted particles as a measure of the temperature - assuming 
that the matter is thermalized - and to look for fluctuations in the event ensemble. Near 
a critical point, the specific heat diverges, causing the thermal fluctuations to vanish. A 
pronounced decrease in the observed temperature fluctuations between otherwise identical 
events would, therefore, indicate an approach to the critical point in the phase diagram. 

The simplest way to identify temperature fluctuations is to look for fluctuations in 
the mean transverse momentum (pt)- Since the increase of Cy near T c is caused by 
the increased fluctuations of the modes associated with the order parameter - in the 
case of the chiral phase transition, the quark condensate (qq) - it should be possible to 
enhance the signal by considering event-by-event fluctuations of observables most sensitive 
to these modes. For hadronic matter, this would be the mean transverse momentum of 
1ow-pt pions, which are thought to partially arise from the decay of local excitations of 
the iso-singlet order parameter, the a-meson mode. 



For practical purposes, two observables have been defined fl2fl : 

1/2 




where is the total squared transverse momentum of all particles emitted in an event, 
N is the event multiplicity, (• • •) denotes the event average, and Ap? the average for the 



inclusive single particle distribution; and |TB 

(N)(AP*/N) 
Ap T 



Fp = ^ v/ ^y iv/ . (6) 



Being a ratio rather than a difference, F p may be less sensitive to collective flow effects 
than $p. An estimate for the expected size of F due to the fluctuations in the cr-field has 
can be obtained in the framework of the linear sigma model pi : 



where £ CT is the correlation length of the fluctuations. 

This prediction has to be modified for two reasons. First of all, the hadronic system 
does not decay into free hadrons right at the critical point, even if it reaches this point 
during its evolution. Instead of observing the large correlation length predicted to occur 
at (T c , fMc) in equilibrium, one would thus expect to observe the reduced correlation length 
associated with the freeze-out parameters (T/,/x/), see Fig. |l|. Furthermore, even under 
conditions of criticality, the correlation length only diverges if the system has enough 
time to develop long-range fluctuations of the order parameter. The theory of dynamical 
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critical phenomena yields an equation determining the change in the inverse correlation 
length m(t) = as the system evolves 



dm 
~dt 



-T(m(t)) (m(t) -m eq (t)) 



(8) 



Here m eq is the inverse correlation length in equilibrium - which vanishes at the critical 
point - and the relaxation rate T(m) ~ m z with some positive exponent z. As critical 
conditions are approached, the relaxation time diverges exhibiting the well-known effect 
of critical "slowing down". 




Figure 1. Schematic QCD phase diagram. The first-order transition line separating the 
quark-gluon plasma (QGP) phase and the hadronic gas (HG) phase terminates in a critical 
point (T c , fi c ), where the transition is of second order. Under favorable conditions, the hot 
system created in a heavy ion collision may transiently pass through the critical point, 
before it freezes out into hadrons that no longer interact. 



If the critical point is reached in a nuclear collision, the system rapidly expands across 
it and freezes out soon afterwards. According to (d) the correlation length £(t) remains 
smaller than £ eq as T c is approached from above, but then decreases less rapidly than £ eq . 
The prediction is thus that the effective correlation length at the moment of freeze-out 
is larger than one would expect under equilibrium conditions, reflecting the temporary 



proximity to the critical point 14,15 



Experimental data for (pt) fluctuations have been obtained by NA49 for Pb+Pb col- 
lisions at the CERN-SPS and recently by STAR for Au+Au collisions at RHIC 



\ \[% . While the NA49 data reflect mean px fluctuations in a forward rapidity region, the 
STAR data were taken at midrapidity. NA49 reported values of F = 1.004 ± 0.004 or 
$ p = 0.6 ± 1 MeV after corrections for two-particle (HBT) correlations. The value from 



STAR is $ p « 35 MeV. 
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4. CHEMICAL FLUCTUATIONS 

Fluctuations in the chemical composition of the emitted hadron yields can potentially 
probe the microscopic structure of the emitting matter. An example is the quark number 
susceptibility 

Xq = T^-{q 1 \) 1 (9) 

which measures the response of the net quark density to a change in the quark chemical 
potenial. A different, but related quantity is the chiral susceptibility Xm = Td(qq) /dm, 
which measures the response of the quark condensate to a change in the current quark 
mass. Xq nas an isoscalar and an isovector component 

Xs = r -1 ([A(« t tt) + A((fid)) 2 ), 
Xns = r -1 ([A(M f ii) - A(d f d)] 2 ), (10) 



which have both been determined on the lattice ||18|| . The results show that xs ~ Xns 



indicating that fluctuations in the u- and d-quark densities are uncorrelated. 

As chemical properties are generally carried by particles, chemical fluctuations are 
determined by fluctuations in the corresponding particle numbers, and their changes are 
governed by particle transport processes. A weakly coupled hadronic gas (HG) and a 
perturbative quark-gluon plasma (QGP) differ significantly in this respect, as listed in 
Table |T|. This implies that measurable differences in the fluctuations associated with 



Table 1 

Quantum numbers of electric charge (Q) and baryon number (B) of particles and their 
abundances in a weakly coupled hadron gas (HG) and quark-gluon plasma (QGP). 





101 


particles 


fraction 


HG 
QGP 


e 

e/3,2e/3 


q,q 


~ 2/3 
~ 1/2 







B 




particles 


fraction 


HG 
QGP 


1 

1/3 


N,N,A,... 

q,Q 


< 1 
~ 1/2 



charge and baryon number exist between a hadron gas and a quark-gluon plasma. 

The second important point is that fluctuations of locally conserved quantities, such 
as the net electric charge or the net baryon number (or the net strangeness!) cannot be 
erased by local reactions. They can only be modified due to particle transport over larger 
distances, i.e. diffusion, and may thus be fixated early, if the evolution of the system is 
sufficiently rapid |19",|20"|. Consider the total net electric charge of all particles emitted 
within a rapidity window Ay in a high-energy nuclear collision. The fluctuation of this 
quantity (AQ 2 )& y is proportional to the volume associated with emission into this rapidity 
window, i.e. to Ay itself in a boost invariant scenario. On the other hand, the rate of 
diffusion of net charge into and out of this volume is independent of the size of the rapidity 
interval. Therefore, early fluctuations of the net charge contained in a sufficiently large 
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rapidity interval Ay > Ay min may survive. The questions are, how large is Ay min , and 
how much do the predictions for a HG and a QGP differ? 

We first address the last question. We estimate the fluctuations per unit volume in the 
two different phases in the weak coupling limit. In the hadron gas we have 

(A5 2 ) HG = N B + N S , 

(AQ 2 ) HG « jV>+jW+iW+iV x - + ..., (If) 

where iVj denotes the number of hadrons of species % in the subvolume, and Nb {N B ) 
counts all (anti-)baryons. In the quark-gluon plasma we have 

(AB 2 )qgp = \(N q + N,), 

(AQ 2 ) QGP « ±(N u + N») + ±(N d + N s + N a + N s ), (12) 

with the analogous notation of the quark numbers. Neglecting the contribution from 
strange quarks, the net charge and baryon number fluctuations in the QGP are related as 

(A,£? 2 )qgp _ 2 

(AQ 2 ) QGP 5- 1 ) 
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Figure 2. Mean square fluctuations per unit entropy of baryon number and charge for 
a weakly interacting hadronic gas versus a weakly interacting quark-gluon plasma. The 
change with beam energy is due to the dependence on baryon chemical potential [I9[ . 



Since it is difficult, if not impossible, to measure the size of an observed subvolume 
accurately, it is desirable to normalize the fluctuations to another extensive quantity. 
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One possibility is to choose the total entropy S, which is closely related to the total 
number of particles in the subvolume. For the QGP one can derive the analytical relation 
(// denotes the baryon chemical potential) 



(A^ 2 )qgp 5 / 22 / fi y 
5 QG p 37tt 2 \ 111 UtJ 



+ ... , (14) 



which also determines the charge fluctuations by virtue of ([13]). Quantitative estimates 
for the hadron gas phase are somewhat more complex, requiring corrections for the decay 
of higher resonances. The resulting estimates for the net charge and baryon number 
fluctuations, divided by the entropy, are compared in Fig. |2| with the QGP predictions 
over an energy range spanning from the SPS to the LHC. The HG prediction is always 
significantly higher, and at RHIC energies and beyond the net charge fluctuations provide 
for greater discrimination than the baryon fluctuations. However, we will argue shortly 
that early net baryon number fluctuations have a greater chance of survival in a narrow 
rapidity interval. 




Ay 

Figure 3. We consider net charge and baryon number fluctuations in a cylindrical volume 
V, defined by a rapidity window Ay in the boost invariant Bjorken model. Conservation 
laws dictate that the total charge and baryon number in V can only change due to particle 
transport through the endcaps of V. 



Next we address the question whether the fluctuations generated during a QGP phase 
can survive the final state expansion in the HG phase. For simplicity, we consider the 
boost invariant case of a longitudinal expansion introduced by Bjorken, illustrated in 
Fig. |3|. Considering a rapidity interval Ay = 1, we see that the left and right boundary 
separate from another with the relative velocity Av ~ 0.76c. The typical longitudinal 
velocity component of a baryon in the hadronic phase is 

v z = (8T/ttM) 1/2 « 0.32c. (15) 

implying that most baryons will not be able to move out of the rapidity interval, in which 
they were contained at the moment of hadronization. Thus, we would expect that the 
baryon number fluctuations generated in a rapidity interval Ay > 1 during the QGP 
phase will remain frozen during the HG phase. A more detailed analysis of the transport 



8 



of baryons through the left- and right-hand surfaces of a cylindrical volume shows that 
baryon number fluctuations from the QGP phase can survive if Ay > v z ||19|| . 

The time evolution of fluctuations can also be studied in the framework of transport 
theory [^l]]. Consider the distribution of particles in rapidity space as a function of proper 
time r: 

n(y,r) = n C(l (y) + f(y,r), (16) 

where f(y, r) describes the deviation from equilibrium. The relaxation toward equilibrium 
is described by the Langevin equation 

df 1 d 2 f 

7T = / + 7^4 + C, (17) 

or r cq dy l 

where r eq is the relaxation time, 7 is the coefficient describing diffusion of particles in 
rapidity, and£ is the noise term. For a locally conserved particle density, r^ 1 = 0, and the 
decay of a fluctuation f(y) is solely governed by diffusion. Making use of the equivalent 
Fokker-Planck equation for the Fourier components of the rapidity fluctuations f(k, r) one 
can show that the Gaussian width af.{r) of the component f(k,r) satisfies the equation 

£0* = -2 7 k 2 (a 2 k - Xk ) (18) 

where Xk = / dr\e lkri (f (y + rj) f (y)) . There is no relaxation of the fluctuations toward equi- 
librium for the k = component due to the conservation law; higher Fourier components 
decay exponentially wih the rate 2k 2, y. 

The diffusion constant 7 is determined by the rapidity transfer in particle collisions 
and the average time between collisions: 7 = (5y co n) 2 /(2r co n). With some additional 
simplifying assumptions one then finds that the mean square fluctuation (AiV 2 ) in a 
rapidity interval Ay decays according to 

(AN 2 ) = Ay ( X o + K 2 - X o)G(Ay diS /Ay)) , (19) 

where G(x) is a universal function. For large times, (AN 2 ) decays slowly toward its 
equilibrium value. Ay 2 m = 2/7(r)(ir describes the minimal rapidity interval for which 
initial state fluctuations survive. Shuryak and Stephanov obtained the estimates [^TJ 

Ayg « 2.2, AygJ « 0.9 (20) 

for the rapidity intervals controlling charge and baryon number fluctuations, respectively. 
However, these estimates are based on collision rates deduced from a hadronic cascade 
model (RQMD) that does not include a deconfined phase, thus possibly overestimating 
the size of the required rapidity windows. 

5. EXPERIMENTAL CONSIDERATIONS AND RESULTS 

For the net charge fluctuations, an elegant way to eliminate the size of the observed 
volume is to consider the fluctuations in the quantity R = N + /N_. These can be related 



to the net charge fluctuations by means of the relation |20 

D = (N ch ) (AR 2 ) = 4(AQ 2 )/(N ch ). (21) 
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It is useful to correct the quantity D for two trivial effects by defining a variable D as |22[ 

D = C^CyD, (22) 

with = (iV + ) 2 /(iV_) 2 accounting for the net average charge in the rapidity interval 
(due to the baryon excess) and C y = 1 — (N ch )^ y /N^^ correcting for the effect of global 
charge conservation. For an uncorrelated thermal pion gas one expects D = 4. Decays 
of higher meson resonances (such as LO,p,rj) reduce the value of D by about 30%. For a 
weakly interacting QGP one anticipates D pa 1. 

Preliminary results for D have been reported by the NA49 collaboration for Pb+Pb 
collisions at 40, 80, and 158 GeV/u beam energy |T(|. A value slightly in excess of D = 4 



was found for all three beam energies and independent of the selected pseudorapidity 
interval 0.3 < An < 3.5. The STAR collaboration has reported first results for net charge 
fluctuations in Au+Au collisions at a/s^vv = 130 GeV in the pseudorapidity window 
| At/ | < 0.7. The value obtained by STAR is D pa 3, as expected for a nearly baryon-free 



hadron resonance gas [17 



6. SUMMARY AND OUTLOOK 

Fluctuations are sensitive probes of the microscopic structure of dense matter. Local 
temperature fluctuations can indicate the presence of a second-order phase transition or 
critical point, while net charge and baryon number fluctuations provide information about 
the particle modes that govern the charge and baryon number transport. In all cases, the 
rapid expansion of the fireball in its last stages is essential for the survival of fluctuations 
established at early times. 

In addition to thermal fluctuations, the fireball may exhibit fluctuations due to the 



initial state ||23[| . These could be probed independently in p+p or p+A collisions. For 
example, such fluctuations may be probes of the quasiclassical coherent glue fields in a 
fast-moving nucleus |23]. Similarly, large nonstatistical fluctuations, such as disoriented 
chiral condensates, can signal the presence of unstable collective modes at some time 
during the course of the heavy ion collision. 

Preliminary data from the NA49 and STAR experiments indicate the presence of both, 
nontrivial pt fluctuations near central rapidity and net charge fluctuations similar to those 
of a pion or hadron resonance gas. The measurement of net baryon number fluctuations 
remains a formidable experimental challenge, since they require the detection of neutral 
as well as charged baryons. However, the case may not be entirely hopeless, as plans by 
the PHENIX collaboration to measure antineutrons with high efficiency demonstrate p5| . 



Fluctuations are a rich topic with many theoretical challenges and opportunities. Lat- 
tice gauge theory can make quantitative predictions for fluctuations near thermal equi- 
librium. The wide array of possible fluctuation observables remains largely unexplored. 
Microscopic models can simulate final state effects on fluctuations and provide useful 
guidance about which observables contain valuable information. There is ample room for 
the exploration of novel analysis strategies, as our picture of the space-time evolution of 
nuclear collisions at RHIC energy comes into focus. One especially promising class of 
observables are the balance functions, which can provide a more differential measure of 
quantum number fluctuations in phase space p6| . The subject also promises a wealth 
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of experimental data, as several RHIC detectors, foremost PHOBOS and STAR, are well 
equipped to explore event-by-event fluctuations in a wide range of observables. 
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